An exact general remeshing scheme applied to physically conservative voxelization

We present an exact general remeshing scheme to compute analytic integrals of polynomial functions over the intersections between convex polyhedral cells of old and new meshes. In physics applications this allows one to ensure global mass, momentum, and energy conservation while applying higher-order polynomial interpolation. We elaborate on applications of our algorithm arising in the analysis of cosmological N-body data, computer graphics, and continuum mechanics problems. We focus on the particular case of remeshing tetrahedral cells onto a Cartesian grid such that the volume integral of the polynomial density function given on the input mesh is guaranteed to equal the corresponding integral over the output mesh. We refer to this as "physically conservative voxelization". At the core of our method is an algorithm for intersecting two convex polyhedra by successively clipping one against the faces of the other. This algorithm is an implementation of the ideas presented abstractly by Sugihara (1994), who suggests using the planar graph representations of convex polyhedra to ensure topological consistency of the output. This makes our implementation robust to geometric degeneracy in the input. We employ a simplicial decomposition to calculate moment integrals up to quadratic order over the resulting intersection domain. We also address practical issues arising in a software implementation, including numerical stability in geometric calculations, management of cancellation errors, and extension to two dimensions. In a comparison to recent work, we show substantial performance gains. We provide a C implementation intended to be a fast, accurate, and robust tool for geometric calculations on polyhedral mesh elements.Comment: Code implementation available at https://github.com/devonmpowell/r3

Voids in cosmological simulations over cosmic time

We study evolution of voids in cosmological simulations using a new method for tracing voids over cosmic time. The method is based on tracking watershed basins (contiguous regions around density minima) of well developed voids at low redshift, on a regular grid of density field. It enables us to construct a robust and continuous mapping between voids at different redshifts, from initial conditions to the present time. We discuss how the new approach eliminates strong spurious effects of numerical origin when voids evolution is traced by matching voids between successive snapshots (by analogy to halo merger trees). We apply the new method to a cosmological simulation of a standard LambdaCDM cosmological model and study evolution of basic properties of typical voids (with effective radii between 6Mpc/h and 20Mpc/h at redshift z=0) such as volumes, shapes, matter density distributions and relative alignments. The final voids at low redshifts appear to retain a significant part of the configuration acquired in initial conditions. Shapes of voids evolve in a collective way which barely modifies the overall distribution of the axial ratios. The evolution appears to have a weak impact on mutual alignments of voids implying that the present state is in large part set up by the primordial density field. We present evolution of dark matter density profiles computed on iso-density surfaces which comply with the actual shapes of voids. Unlike spherical density profiles, this approach enables us to demonstrate development of theoretically predicted bucket-like shape of the final density profiles indicating a wide flat core and a sharp transition to high-density void walls.Comment: 13 pages, 13 figures; accepted for publication in MNRA

Dynamical properties of z ∼4.5 dusty star-forming galaxies and their connection with local early-type galaxies

There is a large consensus that gas in high-z galaxies is highly turbulent, because of a combination of stellar feedback processes and gravitational instabilities driven by mergers and gas accretion. In this paper, we present the analysis of a sample of five Dusty Star Forming Galaxies (DSFGs) at 4 ≲ z ≲ 5. Taking advantage of the magnifying power of strong gravitational lensing, we quantified their kinematic and dynamical properties from ALMA observations of their [C ii] emission line. We combined the dynamical measurements obtained for these galaxies with those obtained from previous studies to build the largest sample of z ∼ 4.5 galaxies with high-quality data and sub-kpc spatial resolutions, so far. We found that all galaxies in the sample are dynamically cold, with rotation-to-random motion ratios, V/σ, between 7 and 15. The relation between their velocity dispersions and their star formation rates indicates that stellar feedback is sufficient to sustain the turbulence within these galaxies and no further mechanisms are needed. In addition, we performed a rotation curve decomposition to infer the relative contribution of the baryonic (gas, stars) and dark matter components to the total gravitational potentials. This analysis allowed us to compare the structural properties of the studied DSFGs with those of their descendants, the local early type galaxies. In particular, we found that five out of six galaxies of the sample show the dynamical signature of a bulge, indicating that the spheroidal component is already in place at z ∼ 4.5

A novel approach to visibility-space modelling of interferometric gravitational lens observations at high angular resolution

We present a new gravitational lens modelling technique designed to model high-resolution interferometric observations with large numbers of visibilities without the need to pre-average the data in time or frequency. We demonstrate the accuracy of the method using validation tests on mock observations. Using small data sets with $\sim 10^3$ visibilities, we first compare our approach with the more traditional direct Fourier transform (DFT) implementation and direct linear solver. Our tests indicate that our source inversion is indistinguishable from that of the DFT. Our method also infers lens parameters to within 1 to 2 per cent of both the ground truth and DFT, given sufficiently high signal-to-noise ratio (SNR). When the SNR is as low as 5, both approaches lead to errors of several tens of per cent in the lens parameters and a severely disrupted source structure, indicating that this is an issue related to the data quality rather than the modelling technique of choice. We then analyze a large data set with $\sim 10^8$ visibilities and a SNR matching real global Very Long Baseline Interferometry observations of the gravitational lens system MG J0751+2716. The size of the data is such that it cannot be modelled with traditional implementations. Using our novel technique, we find that we can infer the lens parameters and the source brightness distribution, respectively, with an RMS error of 0.25 and 0.97 per cent relative to the ground truth.Comment: Submitted to MNRA

A lensed radio jet at milli-arcsecond resolution II: Constraints on fuzzy dark matter from an extended gravitational arc

Using a single gravitational lens system observed at $\lesssim5$ milli-arcsecond resolution with very long baseline interferometry (VLBI), we place a lower bound on the mass of the fuzzy dark matter (FDM) particle, ruling out $m_\chi \leq 4.4\times10^{-21}~\mathrm{eV}$ with a 20:1 posterior odds ratio relative to a smooth lens model. We generalize our result to non-scalar and multiple-field models, such as vector FDM, with $m_{\chi,\mathrm{vec}} > 1.4 \times 10^{-21}~\mathrm{eV}$. Due to the extended source structure and high angular resolution of the observation, our analysis is directly sensitive to the presence of granule structures in the main dark matter halo of the lens, which is the most generic prediction of FDM theories. A model based on well-understood physics of ultra-light dark matter fields in a gravitational potential well makes our result robust to a wide range of assumed dark matter fractions and velocity dispersions in the lens galaxy. Our result is competitive with other lower bounds on $m_\chi$ from past analyses, which rely on intermediate modelling of structure formation and/or baryonic effects. Higher resolution observations taken at 10 to 100 GHz could improve our constraints by up to 2 orders of magnitude in the future.Comment: 5 pages, 2 figures. Accepted in MNRAS Letter

Differential Response to Soil Salinity in Endangered Key Tree Cactus: Implications for Survival in a Changing Climate

Understanding reasons for biodiversity loss is essential for developing conservation and management strategies and is becoming increasingly urgent with climate change. Growing at elevations <1.4 m in the Florida Keys, USA, the endangered Key tree cactus (Pilosocereus robinii) experienced 84 percent loss of total stems from 1994 to 2007. The most severe losses of 99 and 88 percent stems occurred in the largest populations in the Lower Keys, where nine storms with high wind velocities and storm surges, occurred during this period. In contrast, three populations had substantial stem proliferation. To evaluate possible mortality factors related to changes in climate or forest structure, we examined habitat variables: soil salinity, elevation, canopy cover, and habitat structure near 16 dying or dead and 18 living plants growing in the Lower Keys. Soil salinity and elevation were the preliminary factors that discriminated live and dead plants. Soil salinity was 1.5 times greater, but elevation was 12 cm higher near dead plants than near live plants. However, distribution-wide stem loss was not significantly related to salinity or elevation. Controlled salinity trials indicated that salt tolerance to levels above 40 mM NaCl was related to maternal origin. Salt sensitive plants from the Lower Keys had less stem growth, lower root:shoot ratios, lower potassium: sodium ratios and lower recovery rate, but higher δ 13C than a salt tolerant lineage of unknown origin. Unraveling the genetic structure of salt tolerant and salt sensitive lineages in the Florida Keys will require further genetic tests. Worldwide rare species restricted to fragmented, low-elevation island habitats, with little or no connection to higher ground will face challenges from climate change-related factors. These great conservation challenges will require traditional conservation actions and possibly managed relocation that must be informed by studies such as these

TRY plant trait database – enhanced coverage and open access

Plant traits - the morphological, anatomical, physiological, biochemical and phenological characteristics of plants - determine how plants respond to environmental factors, affect other trophic levels, and influence ecosystem properties and their benefits and detriments to people. Plant trait data thus represent the basis for a vast area of research spanning from evolutionary biology, community and functional ecology, to biodiversity conservation, ecosystem and landscape management, restoration, biogeography and earth system modelling. Since its foundation in 2007, the TRY database of plant traits has grown continuously. It now provides unprecedented data coverage under an open access data policy and is the main plant trait database used by the research community worldwide. Increasingly, the TRY database also supports new frontiers of trait‐based plant research, including the identification of data gaps and the subsequent mobilization or measurement of new data. To support this development, in this article we evaluate the extent of the trait data compiled in TRY and analyse emerging patterns of data coverage and representativeness. Best species coverage is achieved for categorical traits - almost complete coverage for ‘plant growth form’. However, most traits relevant for ecology and vegetation modelling are characterized by continuous intraspecific variation and trait–environmental relationships. These traits have to be measured on individual plants in their respective environment. Despite unprecedented data coverage, we observe a humbling lack of completeness and representativeness of these continuous traits in many aspects. We, therefore, conclude that reducing data gaps and biases in the TRY database remains a key challenge and requires a coordinated approach to data mobilization and trait measurements. This can only be achieved in collaboration with other initiatives

Outcomes from elective colorectal cancer surgery during the SARS-CoV-2 pandemic

This study aimed to describe the change in surgical practice and the impact of SARS-CoV-2 on mortality after surgical resection of colorectal cancer during the initial phases of the SARS-CoV-2 pandemic